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Abstract. Tunneling measurements have been carried out on heavily underdoped, slightly over- 
doped and partially Ni-substituted Bi2Sr2CaCu208+a; (Bi2212) single crystals by using a break- 
junction technique. We find that in-plane tunneling spectra below Tc are the combination of in- 
coherent part from the pseudogap and coherent quasiparticle peaks. There is a clear correlation 
between the magnitude of the pseudogap and the magnitude of the superconducting gap in Bi2212. 
The analysis of the data suggests that the tunneling pseudogap in Bi2212 is predominantly a charge- 
density-wave gap on dynamical charge stripes. The tunneling characteristics corresponding to the 
quasiparticle peaks are in excellent agreement with theoretical predictions made for a quasi-one 
dimensional topological-excitation liquid. In addition, the analysis of data measured by different 
techniques shows that the phase coherence along the c-axis is established at Tc due to spin fluctu- 
ations in local antiferromagnetic domains of Cu02 planes. 



I. INTRODUCTION 

Soon after the discovery of superconductivity in 
cuprates , it became clear that the concept of the Fermi 
liquid is not apphcable to cuprates: the normal state 
properties of cuprates are markedly different from those 
of conventional metals (2). The pseudogap (PG) which 
appears in electronic excitation spectra of cuprates above 
Tc, is one of the main features of high- Tc superconductors 
(SCs). There is a consensus on doping dependence of the 
PG in hole-doped cuprates: the magnitude of the PG de- 
creases with increase in hole concentration [^|[^]. How- 
ever, there is a clear discrepancy between the phase dia- 
grams inferred from transport measurements, on the one 
hand, and from tunneling measurements, on the other 
hand: transport measurements |^,^,^ show that, in the 
overdoped region, the PG is absent above Tc, at the same 
time, in tunnehng measurements 

the PG IS 

observed well above Tc. 

The SC characteristics in cuprates have different dop- 
ing dependences: the Tc value has approximately the 
parabolic dependence on hole concentration, p, with the 
maximum near p = 0.16 Q, whereas the SC condensa- 
tion energy has the maximum in the overdoped region 
nearp = 0.19 |Jlltl. 

Recent intrinsic c-axis tunneling data obtained in 
Bi2Sr2CaCu208+x (Bi2212) mesas show that the 
pseudogap (PG) is the normal-state gap, and the PG and 
the superconducting gap (SG) coexist below Tc. High- 
resolution angle-resolved photoemission (ARPES) data 
obtained in Bi2212 at momentum near the (0, tt) 
also show that the quasiparticle (QP) peak and the gap 
have different origins, and they coexist below Tc. Thus, 
the PG in Bi2212 arises from either charge-density waves 
(CDW) or local antiferromagnetic (AF) correlations [or 
spin-density waves (SDW)], or from their combination 
^ , P^JT5l|l^ . It was proposed that, in order to develop 
further understanding of ARPES spectra, it is necessary 



to separate the coherent QP peak from the gap in ARPES 



spectra |11 



Tunneling spectroscopy is an unique probe of SC state 
in that it can, in principle, reveal the QP excitation den- 
sity of states (DOS) directly with high energy resolution. 
In this paper we present tunneling measurements per- 
formed on heavily underdoped, slightly overdoped and 
partially Ni-substituted Bi2212 single crystals by using a 
break-junction technique. We find that in-plane tunnel- 
ing spectra below Tc are the combination of incoherent 
part from the PG and coherent QP peaks. There is a 
clear correlation between the magnitude of the PG and 
the magnitude of the SG in Bi2212. Analysis of the data 
suggests that the tunneling PG in Bi2212 is predomi- 
nantly a CDW gap on dynamical charge stripes. The tun- 
neling characteristics corresponding to the QP peaks are 
in excellent agreement with theoretical predictions made 
for a quasi-one dimensional (ID) topological-excitation 
liquid. In addition, analysis of data measured by dif- 
ferent techniques shows that the phase coherence along 
the c-axis is established at Tc due to spin fluctuations 
in local AF domains of Cu02 planes. In the framework 
of the quasi-lD topological-excitation-liquid scenario in 
Bi2212, the discrepancy between the phase diagrams in 
the overdoped region, inferred from transport and tunnel- 
ing measurements, becomes obvious: they measure two 
different PGs. Other data obtained in cuprates can be 
naturally understood in the framework of such a scenario 
for the SC in Bi2212. To our knowledge, the observation 
of a quasi-lD topological-excitation liquid in Bi2212 is 
presented in the literature for the first time. 



II. THEORY AND MEASURED DATA 

By performing tunneling measurements one can obtain 
the I{V) and dIldV{V) characteristics. The dIldV{V) 
tunneling characteristic measured in a SC-insulator- 
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normal metal (SIN) junction corresponds directly to the 
DOS of QP excitations In this Section, first, we con- 
sider a theoretical I{V) tunneling characteristic in a SIN 
junction and a measured I{V) curve in a SC-insulator-SC 
(SIS) junction of Bi2212. 

Figure 1(a) shows a theoretical I{V) characteristic in a 
SIN junction (Fig. 6 in Ref. p^). In the tunneling legime, 
it is expected that the I{V) curve at high positive (low 
negative) bias, depending on the normal resistance of 
junction, lies somewhat below (above) the normal-state 
curve [the dash line in Fig. 1(a)] [Q. In conventional 
SCs, this prediction is verified by tunneling experiments 
However, we find that, in cuprates, the prediction is 
violated. Figure 1(b) shows the I{V) curve measured in 
a slightly underdoped Bi2212 single crystal with Tc = 83 
K (Fig.l in Ref. In Fig. 1(b), one can see that the 

I{V) curve at high positive (low negative) bias passes 
not below (above) the straight dash line but far above 
(below) the line. This fact cannot be explained by the 
d-wave symmetry of the order parameter. To our knowl- 
edge, this question has been never raised in the literature 
before. This finding is the main motivation of the present 
work. 

Before we discuss our tunneling data, it is necessary to 
consider theoretical characteristics of a topological soli- 
ton and a bound state of two solitons. A topological 
soliton is an extremely stable nonlinear excitation which 
can be moving or entirely static |^,^,|2|,|3|j2|] . The 
topological soliton is often called a kink. A kink is a do- 
main wall between two degenerate ground states. Soli- 
tons and kinks have particlelike properties. Usually, the 
kinks appear in half-filled (spin-) Peierls systems which 
are characterized by either a CDW or SDW instability 
@|l]j22[|2|,|2|,||,|6|,|7). The CDW and SDW instabili- 
ties in (spin-) Peierls systems are the consequence of lat- 
tice distortion. For example, in polyacetylene, the soli- 
tons are moving domain walls between the two degener- 
ate dimerized phases Solitons (kinks) may carry a 
charge (which may be fractional) or spin of 1/2 or both. 
Figure 2(a) shows the characteristics of a kink. Since, in 
Fig. 2, we consider the general solutions, the gap shown 
schematically in Fig. 2 is either a CDW or SDW gap. So, 
the kinks occupy the midgap states. A kink (soliton) can 
form a bound state with another kink (soliton) , which can 
be moving or stationary. Figure 2(b) shows the charac- 
teristics of a bound state of two solitons. 

To our knowledge, a soliton SC was for the first time 
considered in Ref. |^8| in order to explain the SC in 
quasi-lD organic conductors. In this model, two electron- 
solitons form a bisoliton having 2e charge and zero spin 
(a singlet state). The QPs surrounded by deformation 
move coherently as a unique entity along stacks of organic 
molecules without resistance. Later, Davydov |^,^ ap- 
plied the model of the bisoliton SC to cuprates. Re- 
cently, other authors proposed that the charge carriers 
in the normal state of cuprates reside on domain walls 
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FIG. 1. (a) Theoretical /(V^) tunneling characteristic in a 
SIN junction of a SC with the isotropic energy gap [Q. The 
dash line shows a normal-state curve, (b) Measured /(V^) 
curve in a SIS junction of an underdoped Bi2212 with Tc = 
83 K [^. The dash line which is parallel to the liy] curve 
at high bias is a guide to the eye. The arrows show the offset 
from the dash line. Note that the dash line in the plot (b) is 
not the normal-state curve. In fact, any straight line passing 
through (0,0) will show that the behavior of the /(F) curve 
in the plot (b) deviates from the theory |^^ . 



as a consequence of strong short-range electron-electron 
repulsion [p9|. 
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FIG. 2. (a) Normalized dlldV{V) and liy) characteristics 
of a kink, (b) Normalized dl jdViy) and /(V^) characteris- 
tics of a bound state of two solitons po| , ^ . In both insets: 
the gap is shown schematically in grey; the /(V) curve cor- 
responding to the gap is not shown, and the hight of the 
conductance soliton peaks depends on the density of added 
(removed) electrons [|25|j26|. 



III. EXPERIMENT 

Overdoped Bi2212 single crystals and Bi2212 single 
crystals in which Cu is partially substituted for Ni (Ni- 
Bi2212) were grown using a self-flux method. The 
dimensions of the samples are typically of 3x1x0.1 
mm'^. The chemical composition of the Bi-2:2:1:2 phase 
in the overdoped crystals corresponds to the formula 
Bi2Sri gCaCui gOg+x as measured by energy dispersive 
X-ray fluorescence (EDAX). The crystallographic a, 6, c 
values of the overdoped crystals are 5.41 A, 5.50 A and 
30.81 A, respectively. The chemical composition of the 
Bi-2:2:1:2 phase in the Ni-Bi2212 crystals corresponds to 
the formula Bi2Sri.95Cao.95(CuNi)2.o508+x as measured 
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FIG. 3. SIS dI/dV{V) and I{V) measured at 14 K within 
the same underdoped Bi2212 single crystal with Tc = 51 K. 
The dl /dV{V) in both plots are normalized at -400 mV. In- 
set in the plot (b): the differences {dI/dV)a — (dl /dV)b and 
Ia,norm — Ib.norm- The inset in the plot (a) shows how la and 
7;, were normalized: la is normalized at -400 mV, and lb is 
adjusted to be parallel at high bias to la (such procedure is 
equivalent to the normalization at ±00). 

by EDAX. The Ni content with respect to Cu is approx- 
imately 1.5%. The Tc value was determined by four- 
contact method. The transition width is less than 1 K in 
the overdoped Bi2212 crystals, and a few degrees in the 
Ni-Bi2212 single crystals. The underdoped samples were 
obtained from the overdoped single crystals by annealing 
them in vacuum. 

The results presented here are obtained in three sin- 
gle crystals: one underdoped, one overdoped and one 
Ni-Bi2212, having the Tc values of 51, 88 and 75, respec- 
tively. The hole concentration in each of them estimated 
from empirical relation Tc/Tc.max = 1 — 82.6(p — 0.16)^ 
is equal to 0.085, 0.19 and 0.2, respectively. We use 
Tc,max — 95, 95 and 87 K, respectively. The maximum Tc 
value for the Ni-Bi2212 single crystals is estimated from 
dTc/driNi — -5 K/at.%, where n^i is the Ni content with 
respect to Cu . Initially, the tunneling measurements 
in the Ni-Bi2212 single crystals were motivated by the 
presence of magnetic impurities (Ni) in Bi2212. How- 
ever, we use the data obtained in the Ni-Bi2212 as data 
in the overdoped region (p — 0.2) without emphasizing 
the presence of Ni. 

Experimental details of our break-junction setup can 
be found elsewhere | |3l| . In short, many break-junctions 
were prepared by gluing a sample with epoxy on a flexi- 
ble insulating substrate, and then were broken in the ab- 
plane by bending the substrate with a differential screw 



at low temperature in a He ambient. The electrical con- 
tacts (typically with the resistance of a few Ohms) were 
made by attaching gold wires to a crystal with silver 
paint. The I{V) and dI/dV{V) tunneling characteristics 
were determined by the four-terminal method by using 
a standard lock-in modulation technique. At low (con- 
stant) temperature, in one junction, we usually obtain a 
few tunneling spectra by changing the distance between 
broken parts of a crystal, going back and forth etc., and, 
every time, the tunneling occurs most likely in different 
places. 

In addition to SIS-junction measurements, tunneling 
tests have been carried out in the overdoped Bi2212 single 
crystals by forming SIN junctions. Pt-Ir wires sharpened 
mechanically are used as normal tips. 

The magnitude of the SG can, in fact, be derived di- 
rectly from the tunneling spectrum. However, in the ab- 
sence of a generally accepted model for the gap function 
and the DOS in cuprates, we calculate the gap magnitude 
2A as a half spacing between the coherent QP peaks in 
SIS conductance tunneling characteristics. 

IV. RESULTS 

A. The underdoped region 

Figure 3(a) shows the SIS dI/dV{V) and I{V) curves 
obtained in the underdoped Bi2212 single crystal, which 
look like usual tunneling spectra in Bi2212 |Q. In 
Fig. 3 (a), the Josephson /ci?n product is estimated to 
be 13.4 mV. The gap magnitude, Asc = 64 meV, is in 
good agreement with other tunneling measurements . 
The dI/dV{V) and I{V) curves shown in Fig. 3(b) are 
obtained within the same underdoped single crystal as 
those in Fig. 3(a). In Fig. 3(b), the wide humps resem- 
ble the humps in the dI/dV{V) shown in Fig. 3(a). It is 
suggestive that the spectra in Fig. 3(b) correspond to the 
PG. The SC in Bi2212 is weak in the heavily underdoped 
region {p < 0.1) @,|ll|. For example, in YBa2Cu306+i 
(YBCO), modest magnetic fields suppress the SC signif- 
icantly in the heavily underdoped region [^ . This may 
explain why it is possible to observe separately the PG 
in the heavily underdoped Bi2212 by taking into account 
that tunneling spectroscopy probes the local DOS. The 
absence of the Josephson current in the spectra shown in 
Fig. 3(b) indicates that the humps in the conductance are 
incoherent. The differences between the spectra shown in 
Figs 3(a) and 3(b) are presented in the inset of Fig. 3(b), 
which correspond to a "pure SG". Some parts of the 
dI/dV{V) in the inset of Fig. 3(b) are slightly below zero 
because the spectra shown in Figs 3(a) and 3(b) arc not 
taken under the exact same conditions. The small humps 
in dI/dV{V) shown in the inset of Fig. 3(b) are discussed 
below. The dI/dV{V) and I{V) spectra in the inset of 
Fig. 3(b) resemble the characteristics of a bound state of 
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two solitons, shown in Fig. 2(b). Thus, we find that the 
conductance in Fig. 3(a) consists of two contributions: the 
incoherent humps which correspond to the PG (which is 
a normal-state gap 1^), and the coherent QP peaks from 
condensed sohtonhke excitations. In Fig. 3(b), one can 
see that the PG is anisotropic. 

We discuss now the PG. Quasi-ID topological excita- 
tions reside on quasi- ID "objects". Recently, quasi- ID 
charge stripes separated by 2D AF domains have been 
discovered in some cuprates [^,^ . The presence of soli- 
tonlike excitations on stripes implies that the stripes are 
insulating. In other words, there is a charge gap along 
stripes, which is most likely a CDW gap Then, the tun- 
neling PG is predominantly a CDW gap. The same con- 
clusion can be independently obtained from the analysis 
of the data: The magnitude of PG, Ap^ = 130 meV, 
shown in Fig. 3(b), is too large to be explained by the 
development of local AF correlations, since the value of 
the superexchange energy, J ~ 120 meV [Q, is not large 
enough to fit the data. Moreover, by taking into account 
that the PG increases with decrease in hole concentra- 
tion, at lower doping {p < 0.085), the magnitude of the 
PG will become larger than 130 meV. Thus, one can 
conclude that the PG is too large to be the spin gap due 
to local AF correlations. For the second time, indepen- 
dently, we find that the tunneling PG is most likely a 
CDW gap. 

It is possible that the subgap in the dI/dV{V) shown 
in Fig. 3(b) is due to the stripe excitations which were 
damped at higher bias by tunneling electrons. The SC is 
weak in the underdoped region 



B. The overdoped region 

Figures 4(a) and 4(b) show the SIS dI/dV{V) and 
I{V) curves obtained within the same overdoped Bi2212 
single crystal. In Figs 4(a) and 4(b), the Josephson 
IcRn product is estimated to be 6 and 7.5 mV, respec- 
tively. The dI/dV{V) and I{V) characteristics shown in 
Fig. 4(a) look like usual tunneling spectra in Bi2212 [^ . 
The dI/dV{V) and I{V) curves in Fig. 4(b) resemble the 
spectra shown in the inset of Fig. 3(b), and the character- 
istics of a bound state of two solitons in Fig. 2(b). This 
means that, the contribution from the PG in the tun- 
neling spectra shown in Fig. 4(b) is small in comparison 
with the contribution from the QP peaks, at least, at 
low bias. This is most likely due to the fact that, in 
slightly overdoped cuprates, the SC is the strongest, and 
the "strength" of the PG is weak [||jl|]. At high bias, 
the contribution from the PG will be always predomi- 
nant, even if the PG is weak. 

In order to show that we observe not a SIS-junction 
effect but an intrinsic effect, we performed measurements 
in the overdoped Bi2212 single crystals by SIN junctions. 
The inset of Fig.4(a) shows the SIN dI/dV{V) and I{V) 
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FIG. 4. SIS dI/dV(V) and /(V) measured at 14 K within 
the same overdoped Bi2212 single crystal with Tc — 88 K. 
Inset in the plot (a): SIN dI/dV{V) and I{V) measured at 
9 K in an overdoped Bi2212 with Tc = 87.5 K (same axis 
parameters as main plot). The dash lines in the plot (a) and 
in the inset, which are parallel at high bias to the /(V) curves, 
are guides to the eye. Inset in the plot (b): the PG measured 
at 122 K. 



curves obtained in an overdoped Bi2212. In Fig. 4(a), one 
can see that, basically, there is no difference between the 
I{V) characteristics measured in SIS and SIN junctions 
[see the dash lines in Fig. 4(a) and in the inset of Fig. 4(a)] . 

The inset of Fig. 4(b) shows the PG in the slightly over- 
doped Bi2212 at 122 K, which looks like the typical CDW 
gap m quasi- ID conductor NbScg |||,|||. The value of 
122 K is chosen not by accident. Two independent tun- 
neling studies show that, in slightly overdoped Bi2212, 
the onset of SC occurs at 110-116 K @J|l. Thus, the 
value of 122 K is above the onset of SC in slightly over- 
doped Bi2212. In the inset of Fig. 4(b), one can see that 
the PG is anisotropic. The shape of the PG at low bias 
indicates that it has most likely an anisotropic s-wave 
rather than the d-wave shape. 

We discuss now the tunneling data measured in Ni- 
Bi2212 single crystals which are also overdoped in oxy- 
gen. Figures 5(a) and 5(b) show the SIS dl/dViV) and 
I{V) curves obtained within the same Ni-Bi2212 single 
crystal. In Figs 5(a) and 5(b), the Josephson IcRn prod- 
uct is estimated to be 3.1 and 24.5 mV, respectively. 
The SIS dI/dV{V) and liV) curves shown in Fig.5(a) 
look like usual tunneling spectra in Bi2212 ||3^]. The 
dl /dV{V) and I{V) curves in Fig. 5(b) resemble the spec- 
tra in the inset of Fig. 3(b) and the characteristics of a 
bound state of two solitons shown in Fig. 2(b). The ef- 
fect of the absence of the contribution from the PG in the 
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FIG. 5. SIS dI/dV(V) and I{V) measured at 15 K within 
the same Ni-Bi2212 single crystal with Tc = 75 K. Inset in 
the plot (b): dI/dV{V) and I{V) characteristics of the main 
plot at 70.3 K. The dash lines which are parallel to the I{V) 
curves at high bias are guides to the eye. 



tunneling spectra shown in Fig. 5(b) is even stronger than 
that in Fig. 4(b). It is important to emphasize that the 
spectra measured in underdoped Bi2212, shown in in the 
inset of Fig.3(b), and the spectra measured in overdoped 
Bi2212 and Ni-Bi2212, shown in Figs 4(b) and 5(b), are 
similar. So, the data obtained in underdoped and over- 
doped Bi2212 are consistent with each other. However, 
in the overdoped region, the "strengths" of the PG and 
SG are reversed in comparison with those in the under- 
doped region. In the overdoped region, the SC is the 
strongest, and the PG is weak fl|rT|. In Fig. 5(a), one 
can see that the humps in the conductance measured in 
Ni-doped Bi2212 are weaker than those in Figs 3(a) and 
4(a). This means that Ni destroys the PG (CDW gap), 
so stripes become more metallic. This is true from an- 
other point of view: many of our attempts to measure 
the PG above in Ni-Bi2212 have failed. 

Figure 6 displays temperature dependence of the con- 
ductance shown in Fig. 5(b). In Fig. 6, the QP peaks dis- 
appear below Tc — 75 K. This is due to the fact that 
tunneling spectroscopy probes the local DOS; the Tc is 
the macroscopic characteristic. Thus, in Fig. 6, Tc^iocai — 
70 K. Apparently, the data shown in Fig. 6 are measured 
near an impurity {i.e. Ni). It seems that the small humps 
which appear in the dI/dV{V) shown in Figs 4(b) and 
5(b) at bias twice as large as bias of the QP peaks, relate 
to the QP peaks and not to the PG. In Fig. 6, their tem- 
perature dependence resembles the temperature depen- 
dence of a SG. In the inset of Fig. 3(b), a similar hump 
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FIG. 6. Temperature dependence of the conductance peaks 
shown in Fig. 5(b). The conductance scale corresponds to the 
70.3 K spectrum, the other spectra are offset vertically for 
clarity. 



is present in the conductance at negative bias. These 
humps can be well understood in terms of a nanopteron 
soliton which is discussed below. 

The dI/dV{V) and I{V) characteristics at 70.3 K 
in the Ni-Bi2212 are shown separately in the inset of 
Fig. 5(b). If we subtract the linear contribution in the 
I{V) curve shown in the inset of Fig. 5(b) (see the dash 
line), then the dI/dV{V) and I{V) spectra in the inset 
of Fig. 5(b) resemble the characteristics of the kink shown 
in Fig.2(a). 
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FIG. 7. Low temperature phase diagram of Bi2212 
based on the present data: Asc (squares), Aps (dia- 
monds), and Tc (dots) The solid line corresponds to 



Tc/n ^max 
to the eye. 



l-82.6(p-0.16)^ m. The dash lines are guides 



5 



C. The phase diagram 

Figure 7 depicts the phase diagram of Bi2212 based on 
the present data. The magnitude of the SG is in good 
agreement with other tunnehng measurements | |3^ . In 
Fig. 7, one can see that A^c and Tc do not correlate with 
each other. As mentioned above, the tunneling PG shown 
in Fig. 7 is different from the PG inferred from transport 
measurements 

V. THE MEASURED DATA AND THEORY 

We compare now the tunnehng data with theory. The 
characteristics of a kink and a bound state of two soUtons 
are described by hyperboUc functions @|l]j2|j2|,|2|. 
Since the conductance peaks of a bound state of two soh- 
tons, which are shown in Fig. 2(b), look very similar to 
the conductance peaks not only of high- Tc SCs but also 
of low-Tc SCs (not the background) , we rely here exclu- 
sively on the I(V) characteristics which are conceptually 
different for the two models based on the Fermi liquid, in 
the first model, and on a quasi-lD topological-excitation 
liquid, in the second model. 

A. I{V) characteristics 

Figure 8(a) shows the measured I{V) curve from the 
inset of Fig. 3(b). In Fig. 8, for simplicity, we analyze 
the data only at positive bias. As shown in Fig. 8(a), 
the data from the inset of Fig. 3(b) can be fitted very 
well by the hyperbolic function f{V) — Ax{tanh[{eV - 
2A)/eVb] + tanh[{eV + 2A)/eVo]), where e is the elec- 
tron charge; V is the bias; A is the maximum SC energy 
gap, and A and Vq are the constants. In Fig. 8(a), we 
also present the measured I{V) curve of the PG from 
Fig. 3(b). We find that any tunneling I{V) curve in 
Bi2212 can be resolved into the two components shown 
in Fig. 8(a). Thus, the I{V) characteristics in Bi2212 
(as well as the dI/dV{V) characteristics) consist of two 
contributions: from the quasi-lD topological excitations 
and from the PG which is the normal-state gap The 
"usual" liV) and dI/dV{V) spectra in Bi2212 show the 
presence of both components [see Figs 3(a), 4(a) and 
5(a)]. The absence [see Fig. 3(b)] or weak contribution 
of one component [see Figs 4(b) and 5(b)] in tunneling 
spectra make the appearance of the spectra " unusual" . 

First, we analyze the tunneling data from Ref. [^: 
Figure 8(b) shows the data from Fig. 1(b), the f{V) fit, 
and their difference. Figure 8(c) depicts the two com- 
ponents in the I{V) curve from Fig. 4(b). As shown in 
Fig. 8(d), the contribution from the QP peaks in the SIN 
I{V) from the inset of Fig. 4(a) seems to be weaker than 
that in SIS junctions. As seen in Fig. 8, all plots are sim- 
ilar. To fit the SIN I{V), we use the same f{V) function 




FIG. 8. Measured I{V) curves and the f{V) fit (see text) 
in (a)-(b) underdoped and (c)-(d) overdoped Bi2212: (a) The 
data (diamonds) from the inset of Fig. 3(b); the measured 
I{V) of the PG from Fig.3(b), and the f{V) fit. (b) The 
data from Ref. jl9|, shown in Fig. 1(b); the f{V) fit, and their 
difference, (c) The I{V) from Fig.4(b); the f{V) fit, and their 
difference, (d) The SIN liV) from the inset of Fig.4(a); the 
f{V) fit [2A ^ A in f{V)], and their difference. In the plots 
(b) and (c), the Josephson currents are removed. 

by substituting 2 A for A. In Figs 8(b)-8(d), the am- 
plitude, A, of the f{V) fit can be changed, this only af- 
fects the scale but not the shape of the differences which 
correspond to the PG. The decomposition of the I{V) 
curves in Figs 8(b)-8(d) into the two components ex- 
plains why I{V) tunneling characteristics in cuprates do 
not obey the theoretical predictions based on the Fermi- 
liquid model The I{V) curves in Figs 4(a) and 5(a) 
and in Refs [|l^ and [Q can be resolved into the two 
components in the same manner. The data in Fig. 3 are 
self-consistent. 




FIG. 9. Measured /(V) curves and the f{V) and f„{V) 
fits (see text) in Ni-Bi2212: (a) The I{V) from Fig.5(b); the 
f{V) fit, and their difference, (b) The I{V) from the inset 
of Fig. 5(b) (diamonds) (the linear background is subtracted) 
and the fn{V) fit. In the plot (a), the Josephson current is 
removed. 
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However, we find that the I{V) curve shown in Fig. 5(b) 
consists of three components. Figure 9(a) shows the 
I{V) from Fig. 5(b), the f{V) fit and their difference. In 
Fig. 9 (a), it is clear that the difference consists of two con- 
tributions: the contribution at high bias is predominantly 
from the PG, similar to those in Fig. 8, and the contribu- 
tion at low bias with the peak at 35 mV. By analyzing the 
data we show below that this peak corresponds to elec- 
tron tunneling assisted by spin excitations . In other 
words, this peak corresponds to the magnetic resonance 
peak observed in inelastic neutron scattering experiments 

Figure 9(b) shows the I{V) curve from the inset of 
Fig. 5(b), where the linear contribution [see the dash line 
in the inset of Fig. 5(b)] is subtracted. As mentioned 
above, the I{V) curve in Fig. 9(b) is very similar to the 
I{V) characteristic of a kink, shown in Fig. 2(a). This 
means that, at this temperature, the quasi-lD topologi- 
cal excitations are already decoupled. We use the hyper- 
bolic function fniV) = A\ xtanh{V/Vo) to fit the data in 
Fig. 9(b), which gives a good interpolation. Thus, there 
is good agreement between the data below Tc and above 
(local) Tc- 

In conventional SCs, tunneling conductances measured 
in SIN and SIS junctions contain different information. If 
a SIN conductance corresponds directly to the SC DOS, 
a SIS conductance is the convolution of the DOS with 
itself [Q. In Fig. 8, we use the same function to fit I(y) 
characteristics measured in SIS and SIN junctions. Phys- 
ically, it is incorrect. However, the main point of the fit is 
that the asymptotics of I{V) characteristics correspond- 
ing to QP peaks, measured either in SIS or SIN junctions, 
are "flat" (constant). The "flat" asymptotics in tunneling 
spectra are a hallmark of the presence of one dimension- 
ality in the system. Such asymptotics never appear in 
2D or 3D systems. 

So, we conclude that, in Bi2212, the I{V) character- 
istics corresponding to the QP peaks definitely disagree 
with Blonder-Tinkham-Klapwijk's predictions based on 
the Fermi-liquid model fl^ , and are in good agreement 
with the theoretical predictions made for a quasi- ID 
topological-excitation liquid [p0pl|j2^,p^,p3 . 



B. dI/dV{V) characteristics 

As emphasized above, we do not rely on the dI/dV{V) 
data in order to draw any conclusion about the origin of 
QPs in Bi2212. However, the analysis of the data would 
be not complete, if we did not consider the conductance 
spectra. The coherent QP peaks in the SIS tunneling 
spectra can be fitted very well by the derivative [f{V)]' — 
A2x{{sech[{eV -2A)/eVo])^ + {sech[{eV + 2A)/eVb])^), 
as shown in Fig. 10. The fit is applicable only to the coher- 
ent QP peaks, and not to the humps at high bias, which 
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FIG. 10. Measured dI/dV{V) curves (circles) and the 
[f{V)]' fit (see text): (a) dI/dV{V) from the inset of Fig.3(b); 
(b) dI/dV{V) from Fig.4(b); (c) dI/dV{V) from Fig.5(b); (d) 
SIN dI/dV{V) from the inset of Fig.4(a) (2A ^ A in [f{V)]'). 



correspond to the PG, and neither to the Josephson cur- 
rent. The dI/dV(V) curves in Figs 4(a) and 5(a) and 
in Refs [|l9| and can be fitted in the same manner. 
For example, the conductance in an overdoped Bi2212 
single crystal, shown in Fig. 2 of Ref. ||l9|, clearly is the 
combination of QP peaks and humps. 

In Fig. 10(c), the coherent QP peaks contain also a con- 
tribution from electron tunneling due to spin excitations 
[see Fig. 9(a)]. As seen in Figs 5(b) and 6, the tops of the 
QP peaks have a composite structure which we attribute 
to the presence of the contribution from electron tunnel- 
ing due to spin excitations: This contribution is literally 
developed on top of the QP peaks. In Fig. 10(d), in order 
to fit the SIN dI/dV{V) at low bias, we had to shift up 
the [fiV)]' curve. 

The conductance shown in the inset of Fig. 5(b) can be 
fitted by [fniV)]' = A3x[sechiV/Vo)]^. 



C. Small humps 

The small humps which appear in the dI/dV{V) at 
bias twice larger than bias of the QP peaks, shown in 
Figs 4(b) and 5(b) as well as in Figs 10(b) and 10(c), can 
be understood in terms of a nanopteron soliton |20|] . A 
nanopteron is a bound state resulting from the nonlinear 
interaction between the soliton (or kink) and the peri- 
odic wave (see Fig. 6. 27(c) in Ref. [|2^). Instead of flat 
asymptotics, the characteristics of a nanopteron have os- 
cillations. A bound state of two solitons can also exist in 
resonance with small amplitude linear waves (see Fig.l 
in Ref. [0). Since, in Figs 4(b) and 5(b), the predomi- 
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nant contribution in the conductances at high bias arises 
from the PG, it is impossible to follow the oscillations in 
the conductances at high bias. Fortunately, the oscilla- 
tions can be also seen in the I{V) characteristics shown 
in Figs 8(a) and 8(c). In Fig. 8(a), the oscillations can be 
observed up to 480 mV. In Fig. 8(c), only one hump is 
present in the I{V) characteristic at about 105 mV. The 
next question is to understand what periodic waves in 
Bi2212 can interact with the bound states of two solitons 
and cause these oscillations? Phonons are the primary 
candidate on this role. It is also possible that magnonlike 
excitations (spin 1, charge 0) which cause the appearance 
of the magnetic resonance peak in neutron spectra ||4C| ] 
correspond to these linear waves, which are in resonance 
with bound states of topological solitons. 

It is well known that dI/dV{V) characteristics mea- 
sured in SIS junctions often have subgap structures [see, 
for example, Fig. 4(a)]. These subgap structures can be 
also understood in terms of a nanopteron soliton. They 
most likely have the same origin as the oscillations at 
high bias, i.e. due to the nonlinear interaction between 
a bound state of two solitons and the periodic wave (see 
Fig.4 in Ref. 



VI. STRIPES AND TOPOLOGICAL 
EXCITATIONS 

Here we discuss the stripes in cuprates and stripes ex- 
citations. 

Analysis of the data obtained in inelastic neutron scat- 
tering experiments suggests that the stripes in cuprates 
are of the 2kp type |4^. The chains in YBCO also have 
the 2kp charge modulation, as observed by tunneling 
spectroscopy |43 . The 2kp ordering pattern on a stripe 
is shown schematically in Fig. 11 (a). Recent simulations 
of ARPES data in Bi2212 show that the stripes in Bi2212 
are most likely site-centered [Q . So we consider only the 
2kp site-centered stripes. 

The stripes in cuprates are dynamical |^^. Recently, 
the dynamics of insulating charge stripes was studied 
by Zaanen and co-workers [^ , ^5p^ ]. Figures 11(b) and 
11(c) show a soliton and a kink on 2kp stripes, respec- 
tively. In the 2kp charge modulation along stripes, a kink 
(soliton) carries an e charge, consequently, a bi-soliton 
has 2e charge. Apparently, charge carriers in nanotubes 
are also topological solitons Q . 

The soliton SC scenario, probably, is realized in quasi- 
ID organic metals [^Sj, ladders [Q, and on CuO chains 
in YBCO which have a CDW ground state @. The lat- 
ter fact can be naturally understood since the stripes in 
Cu02 planes have also the CDW ground state. Recently, 
midgap solitons have been observed on CuO chains in 
YBCO by tunneling spectroscopy j5^. The spectrum 
averaged along a CuO chain shows that there is a weak 
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FIG. 11. 2kF half-filled stripes: (a) reference state; (b) soli- 
ton (S) and antisoliton, and (c) kink (K) and antikink. The 
full/empty circle denotes the presence/absence of the hole. 



bound state of solitons inside the CDW gap. The mag- 
nitude of the induced SG of the bound state of solitons 
is about 6 mcV. 

Charged solitons repel each other. In conventional 
SCs, two electrons which form a Cooper pair also re- 
pel each other. The occurrence of an attractive potential 
between electrons is central to the SC state. So, charged 
solitons will couple with each other if there is an attrac- 
tive potential between them. 

The next question is to understand how is the long- 
range phase coherence established? In Cu02 planes, 
the phase coherence can occur, for example, due to the 
Josephson coupling between SC stripes. Indeed, the so- 
called Yamada plot suggests that Tc increases if and only 
if stripes move closer together Since the SC in Cu02 
planes is fully 2D because the coupled excitations reside 
on stripes, the phase coherence along the c-axis is most 
likely established due to a different mechanism. 



VII. C-AXIS PHASE COHERENCE AND SPIN 
FLUCTUATIONS IN CUPRATES 

There is a consensus that the SC in cuprates is two- 
dimensional (2D). It is widely believed that the long- 
range phase coherence occurs at Tc due to the Josephson 
couphng between SC Cu02 (bi-, tri-, ...)layers. Recent 
measurements of the in-plane (pab) and out-of-plane (pc) 
resistivities in Tl2Ba2CaCu208 as a function of applied 
pressure show that Pc{T) shifts smoothly down with in- 
crease of pressure, however, Tc first increases and then 
decreases js^ . This result can not be explained by the 
interlayer Josephson-coupling mechanism. The authors 
conclude ||5^: "Any model that associates high- Tc with 
the interplane Josephson coupling should therefore be re- 
visited." 

Here we analyze data obtained in Andreev reflec- 
tion, inelastic neutron scattering (INS), microwave, muon 
spin relaxation (/iSR), tunneling and resistivity mea- 
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surements performed on different cuprates, mainly, on 
YBCO, Bi2212 and La2-:,Sr:,Cu04+y (LSCO). Analy- 
sis of the data shows that the long-range phase coher- 
ence in the cuprates intimately relates to AF interactions 
along the c axis. We also analyze data measured in heavy 
fermions UPt3, UPd2Al3 and Celrlns, and in some lay- 
ered non-SC compounds with ferromagnetic (FM) corre- 
lations. 

A. Introduction 

SC and magnetism were earlier considered as mutually 
exclusive phenomena. Recent research revealed a rich va- 
riety of extraordinary SC and magnetic states and phe- 
nomena in novel materials that are due to the interaction 
between SC and magnetism Q. The coexistence of SC 
and long-range AF order was first discovered in /JMoeSes 
{R = Gd, Tb and Er), i?Rh4B4 (R = Nd, Sm and Tm), 
and i?M06S8 {R = Gd, Tb, Dy and Er) [||. Later, coex- 
istence of SC and AF order was found in U-based heavy 
fermions (UPtg, URu2Si2, UNi2Al3, UPd2Al3, UgCo, 
and UgFe), in heavy fermions -RRh2Si2 {R = La and 
Y), Cri^x^Gx, CeRu2, in borocarbides i?Ni2B2C {R — 
Tm, Er, Ho, Dy), in organic S Cs [l^JS^ and in the new 
heavy fermion CeRhi_2;Ira;In5 ]55|]. In CeRho.5Iro.5In5, 
the bulk SC coexists microscopically with small-moment 
magnetism (< O.l/is) |^^. In all other heavy fermions, 
there are strong AF correlations present in the SC state 
p6| , p4|j55| . Another class of materials in which SC and 
AF coexist are cuprates such as YBCO and LSCO com- 
pounds iQ . Coexistence of SC and FM order is found in 
the heavy fermion UGe2 |5^], and in Ru-based materials, 
for example, in RuSr2GdCu2 08 | |57| |. 

In SC heavy-fermion systems, spin fluctuation 
(electron-electron interactions) are believed to mediate 
the electron pairing that leads to SC |5^. For the heavy 
fermions Ceing, CePd2Si2 ||, UPd2Al3 @ and UGe2 
, there is indirect evidence for spin-fluctuation mech- 
anism of SC. This intimate relationship between the SC 
and magnetism also appears to be central to the SC 
cuprates js^ , which inherited the magnetic properties 
from their parent compounds, AF Mott insulators. Many 
theoretical studies suggest that the SC in cuprates is me- 
diated via the exchange of AF spin fluctuations |6|] . 

In cuprates, there are two energy scales the pair- 
ing energy scale, Aj, (see A^c in Fig. 7), and the phase- 
coherence scale, Ac, observed experimentally p2| , |63t . 
The two energy scales have difl^erent dependences on hole 
concentration, p, in Cu02 planes: Ap increases linearly 
with decrease in hole concentration, whereas Ac has ap- 
proximately the parabolic dependence on p and scales 
with Tc as 2Ac ~ bAkeT^ [||. 




T/T 



FIG. 12. Temperature dependences of in-plane (dots) 
B.nd c-axis tunneling QP peaks (squares) [^] in 
slightly overdoped Bi2212 single crystals, A{T) / A{T„ii„) . 
The BCS temperature dependence is shown by the solid line. 
The dashed lines are guides to the eye. 

B. In-plane and c-axis tunneling in Bi2212 

Figure 12 shows the temperature dependence of in- 
plane tunneling QP peaks, measured in slightly over- 
doped Bi2212 single crystals [@,|j|,|l^. Tunneling mea- 
surements performed on slightly underdoped Bi2212 sin- 
gle crystals show similar temperature dependence [ p2[ . 
So, the temperature dependence of in-plane QP peaks 
shown in Fig. 12 can be considered as typical. In Fig. 12, 
we also present the temperature dependence of c -axis QP 
peaks measured in slightly overdoped Bi2212 mesas 
Measurements performed on micron-size mesas present 
intrinsic properties of the material. In Fig. 12, surpris- 
ingly, the temperature dependences along the ab planes 
and along the c axis are different. It is a clear hall- 
mark of the coexistence of two different SC mechanisms 
in Bi2212: in-plane and along the c axis. 

C. LSCO 

In LSCO, there is evidence that the SC intimately re- 
lates to the establishment of AF order along the c axis. 
Recent /iSR measurements performed on non-SC Eu- 
doped LSCO having different hole concentrations show 
that the SC phase of pure LSCO is replaced in Eu-doped 
LSCO by the second AF phase (see Fig.4 in Rcf. @). 
Thus, the data show that it is possible to switch the en- 
tire hole concentration dependent phase diagram from 
SC to AF. It is a clear hallmark that the SC in LSCO in- 
timately relates to the formation of AF order. We return 
to this important result later. 

We turn now to the analysis of resistivity data mea- 
sured in non-SC 2D layered compounds with AF or FM 
correlations. The data clearly show that, in all these 
layered compounds, the out-of-plane resistivity, pc, has 
drastic changes either at Neel temperature, T/v , or Curie 
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temperature, Tc, whereas the in-plane resistivity, pab, 
passes through T/v or Tc smoothly. 

We start with YBCO. In AF undoped YBCO (x = 
0.35; 0.33; and 0.32) having Tn ~ 80 K, 160 K, and 210 
K, respectively, pc shows a sharp increase, by about 2 or- 
ders of magnitude, upon cooling through T/v [|5j. At the 
same time, pat changes smoothly at T/y. The same effect 
has been observed in LuBCO (x = 0.34) @. So, the 
Neel ordering in undoped YBCO has remarkably differ- 
ent impact on the electron transport within Cu02 planes 
and between them. The authors conclude |^ : "The Neel 
temperature actually corresponds to the establishment of 
AF order along c axis." 

We now discuss the FM layered compounds. The struc- 
ture of the FM compound Bi2Sr3Co20g (BSCoO) is sim- 
ilar to the structure of Bi2212, where C0O2 planes are 
analogous with Cu02 planes in Bi2212 |67|. BSCoO be- 
comes FM at Tc ~ 3.2 K. The resistivity data show that, 
at Tc, there is a cusp in pc but pab changes smoothly. The 
authors conclude that the long-range magnetic order in 
BSCoO develops at Tc along the c axis In layered 
manganite Lai.4Sri.gMn207 (LSMO) which is composed 
of the Mn02 bilayers becomes FM at Tc = 90 K [|8). 
The resistivity data show that, at Tc, there are drastic 
changes in pc (a few orders of magnitude) , but very small 
changes in pab- They conclude that, in LSMO, the long- 
range magnetic order develops at Tc = 90 K along the c 
axis Q. 

Neutron-scattering measurements performed on the 
heavy fermion URu2Si2 (Tc = 1.2 K) show that the AF 
order develops at T/v = 17.5 K along the c axis So, 
it seems that, in all layered compounds, the long-range 
AF or FM order develops at T/v or Tc along the c axis 
(the in-plane magnetic correlations exist above T/v and 
Tc M)- 

We return now to the analysis of the phase diagram 
of non-SC Eu-doped LSCO, where the SC phase of pure 
LSCO is replaced by the second AF phase Q. The 
conclusion made in the previous paragraph signifies that 
either the main AF phase of Eu-doped LSCO or the sec- 
ond AF phase develops along the c axis. Thus, the SC 
phase of pure LSCO is replaced in Eu-doped LSCO by the 
AF phase which develops along the c axis. Consequently, 
the SC in LSCO intimately relates to the establishment 
of the long-range AF order along the c axis. 

In heavy fermion Celrlns, /iSR measurements discov- 
ered the onset of a small magnetic field (~ 0.4 Gauss) 
which sets exactly at Tc |^^. In YBCO (x = 0.6), recent 
INS measurements identified small magnetic moments di- 
rected along the c axis, which increase in strength at and 
below Tr fzOtl. 



D. YBCO, Bi2212 and LSCO 

Here we compare coherence SC and magnetic charac- 
teristics of YBCO, Bi2212 and LSCO. The comparison 
shows that the magnetic and coherence SC characteris- 
tics have similar temperature dependencies, and, at dif- 
ferent dopings, their magnitudes are proportional to each 
other (and proportional to Tc). Thus, the coherence SC 
and magnetic properties of cuprates intimately relate to 
each other. 

First, we describe the magnetic properties of cuprates. 
The low energy magnetic excitations in LSCO cuprate 
have been extensively studied, and the observed spin 
fluctuations are characterized by wave vector which is 
incommensurate with the lattice [Q. These modulated 
spin fluctuations in LSCO persist in both normal and 
SC states. The spin dynamics in YBCO and Bi2212 
studied by INS exhibit below Tc a sharp commensu- 
rate resonance peak which appears at well defined en- 
ergy Er ||,^,g7|J§,0,^,|9E,|l|,|| . Incommensu- 
rability in YBCO has been also reported |7^, and it is 
consistent with that in LSCO of the same hole doping, 
but, in YBCO, it occurs in the SC state. Now it is clear 
that the incommensurability and the commensurate res- 
onance are inseparable parts of the general features of the 
spin dynamics in YBCO at all doping levels |^3|. Thus, 
there is a clear evidence of coexistence of AF order and 
SC below Tc, at least, in LSCO and YBCO. 

We now compare coherence SC and magnetic charac- 
teristics of the cuprates. Figure 13(a) shows the temper- 
ature dependences of the superfluid density in near op- 
timally doped single crystals of Bi2212 ^ and YBCO 
|B5|, and in an overdoped LSCO (x = 0.2) single crystal 
1 86 1, measured by microwave, /iSR and ac-susceptibility 
techniques, respectively. The superfluid density is pro- 
portional to 1/A^(T), where A(T) is the magnetic pen- 
etration depth. Figure 13(b) shows the temperature 
dependences of Andrcev-reflection gap measured in an 
overdoped Bi2212 thin film |Q, and in overdoped single 
crystals of Bi2201 (88) and LSCO (x = 0.2) It is 

important to emphasize that Andreev reflections are ex- 
clusively sensitive to coherence properties of condensate. 
In Figs 13(a) and 13(b), one can see that there is good 
agreement among temperature dependences of coherence 
SC characteristics of different cuprates. 

Figure 13(c) shows the temperature dependences of the 
peak intensity of the incommensurate elastic scattering in 
LSCO (x = 0) l?^] and the intensity of the commensurate 
resonance peak measured by INS in near optimally doped 
Bi2212 H and YBCO 0. 

In Fig. 13, all temperature dependences of coherence 
SC and magnetic characteristics below Tc exhibit a strik- 
ing similarity. Since all temperature dependences shown 
in Fig. 13 are similar to the temperature dependence of 
c-axis QP peaks in Bi2212, shown in Fig. 12, and dif- 
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ferent from the temperature dependence of in-plane QP 
peaks, also shown in Fig. 12, it is reasonable to assume 
that the coherence SC and magnetic properties of the 
cuprates intimately relate to each other along the c axis. 
Moreover, in YBCO (x = 0.6), recent INS measurements 
found small magnetic moments directed along the c axis, 
which increase in strength at and below Tc [Q. At the 
same time, as one can see in Figs 12 and 13, the in-plane 
mechanism of the SC in the cuprates has no or little re- 
lations to the magnetic interactions, at least, along the c 
axis. 

Figure 14 shows the energy position of the magnetic 
resonance peak, Er, in Bi2212 ||7|,|7|l and in YBCO 
||,^j76[|7|,|7|,^,|^,|l]j8| as a function of doping. The 
parabolic curve corresponds to the coherence energy 
scale, Ac, which is proportional to Tc as 2Ac ~ 5.4fcs7c 
p2| . At different doping levels, Andreev-reflection data 
coincide with Ac |6^ , |6^ , p9[ . In Fig. 14, we present also the 
A.,^ data from Fig. 7, which correspond to the pairing en- 
In Fig. 14, one can see that, at different 
is proportional to Ac as Er ~ 2 Ac. This 



ergy scale |62 
dopings, Er 

correlation suggests that spin excitations are responsi- 
ble for establishing the phase coherence in YBCO and 
Bi2212 since the relation Er ~ 2Ac is in good agreement 



FIG. 13. (a) Temperature dependence of the superfluid 
density in near optimally doped single crystals of Bi2212 (Tc 
= 93 K) |84| and YBCO {Tc = 93 K) and in an overdoped 
LSCO (x = 0.2) single crystal {T = 36 K) |86|. The super- 
fluid density is proportional to 1/A^(r), where A(r) is the 
magnetic penetration depth. Inset: temperature dependence 
of the superfluid density (squires) in heavy fermion Celrlns 
{Tc = 0.4 K) (axis parameters as main plot), (b) Temper- 
ature dependence of Andreev-reflection gap, A{ T) / A(Tmin) . 
in an overdoped Bi2212 thin fllm (Tc — 80 K) |p7|, and in over- 
doped single crystals of Bi2201 (Tc = 29 K)^^ and LSCO 
(x — 0.2) (Tc — 28 K) ji9|. Inset: temperature dependence 
of the Andreev gap (dots) in heavy fermion UPts {T ~ 440 
mK) (axis parameters as main plot), (c) Temperature de- 
pendence of the peak intensity of the incommensurate elastic 
scattering in LSCO (x = 0) {Tc = 42 K) (ti) and the intensity 
of the magnetic resonance peak measured by INS in near op- 
timally doped Bi2212 {T = 91 K) (7| and YBCO {T = 92.5 
K) |77|. The neutron-scattering data are average, the real 
data have the vertical error of the order of ±10% Q^jT?). 
The BCS temperature dependence is shown by the thick solid 
line. The dashed lines are guides to the eye. 



with the theories in which the SC is mediated by spin 
fluctuations |6l| . Recently, it was shown that applied 
magnetic fields suppress the magnetic resonance peak in 
YBCO, indicating that the resonance peak indeed mea- 
sures the long-range phase coherence . The strength 
of coupling between spin excitations and charge carriers 
is sufficient to account for the high Tc value in cuprates 

What is interesting is that all temperature depen- 
dences shown in Fig. 13 are similar to the temperature 
dependence of the superfluid density in heavy fermion 
Celrln5, measured by /iSR and to the tempera- 

ture dependence of the Andreev-reflection gap in heavy 
fermion UPta |Q , which are shown in the insets of Figs 
13(a) and 13(b), respectively. Spin fluctuations are be- 
lieved to mediate the electron pairing in Celrlns |5^] and 
UPta that leads to SC. The magnetic resonance peak 
has not yet been detected in Celrlns or UPta, however, 
the magnetic resonance peak has been observed by INS in 
another heavy fermion UPd2Al3 where spin fluctua- 
tions mediate the SC which coexists with the long-range 
AF order Q . The latter facts also point to the presence 
of the spin-fluctuation coupling mechanism in cuprates. 

The behavior of all temperature dependences shown 
in Fig. 13 can be easily understood in terms of the 
spin-fiuctuation SC mechanism (electron-electron inter- 
actions): crudely speaking, they exhibit the squared BCS 
temperature dependence. 



E. Discussion 

In spite of the unmistakable similarities among the 
magnetic and SC properties of YBCO, Bi2212 and LSCO 
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(and some heavy fermions for which there is an in- 
direct evidence for spin-fluctuation mechanism of SC), 
clearly, there is a difference between magnetic proper- 
ties of LSCO and YBCO. If, in YBCO, = T^om = 
Tine, where Tcom {Tine) is the onset temperature of the 
(in-) commensurate peak(s), in LSCO the situation is 
different. First, the commensurate peak has not been 
detected. Second, in LSCO, mainly, Tc < Tine (Figure 
13(c) shows the case when Te Tine)- So, the magnetic 
and SC properties of LSCO are similar to those of most 
heavy fermions: the commensurate peak has not been 
detected, and Te < Tjy. For example, in UPd2Al3, Te ~ 
2 K and Tn = 14.5 K @,||; in URuaSiz, Te « 1.2 K 
and T/v = 17.5 K and, in FM heavy fermion UGe2, 
Tc < Tc @. In Bi2212, the situation looks more like 
that of YBCO, even if, the incommensurate peaks have 
not yet been detected. 

In fact, the differences between magnetic and SC prop- 
erties of YBCO and LSCO (and among some heavy 
fermions) can be understood in terms of the following 
chain of events: the formation of pairs at Tpeiij- - the SC 
order parameter couples to the magnetic in-plane order 
parameter - the appearance of AF interactions along the 
c axis at - the appearance of the long-range SC phase 
coherence. If there is no pairs, the SC is absent, even if, 
the AF order is established (as in the Eu-doped LSCO 
p^). If Tpair < Tni, then, this reminds the situation in 
LSCO where Te < Tine- If Tpair > T„i, then, this re- 
minds the case of YBCO where Te — Team = Tine, and, 
consequently, Te = r,„. Unfortunately, this simple pic- 
ture cannot explain the presence of the commensurate 
peak in INS spectra but, at least, shows that the differ- 
ences between the magnetic and SC properties of LSCO 
and YBCO can be understood in the framework of one 
picture with different initial parameters. 



F. Summary 

Analysis of the data obtained by different techniques 
in YBCO, Bi2212 and LSCO shows that the long-range 
phase coherence intimately relates to AF interactions 
along the c axis. Apparently, in cuprates, the magnetic 
and SC order parameters are coupled to each other, and 
the phase-coherence scale, Ac, has the magnetic origin 
(see Fig. 14). 

There are common features in the SC state of the heavy 
fermions Celrln5, UPta and UPd2Al3, on the one hand, 
and the cuprates, on the other hand. It is possible that, in 
all heavy-fermion and organic SCs, the long-range phase 
coherence is established due to spin fluctuations. The 
pairing mechanism may be different (as in Bi2212). 




0.08 0.16 0.24 0.32 

hole concentration, p 

FIG. 14. Phase diagram: the energy position of the mag- 
netic resonance peak , Er , in Bi2212 (squares) [ {7^ , [73| and 
in YBCO (dots) 0,^|j7|,0j7|^,||,|l|]ii and A^e (dia- 
monds) from Fig. 7. The parabolic solid line corresponds to 
the coherence energy scale, Ac which is proportional to 
Tc as 2Ac — dAksTc. When not available, p has been cal- 
culated from the relation Tc/Tc,max = 1 - 82.6(p - 0.16)^ §], 
and we use Tr 



= 93 K and T^, 
YBCO and Bi2212, respectively, 
the eye. 



= 95 K in the case of 
The dash line is a guide to 



VIII. TUNNELING ASSISTED BY SPIN 
EXCITATIONS IN BI2212 

In Fig. 14, one can see that at p = 0.2, Er ~ 2Ase- 
Therefore, the tunneling in a SIS junction of Bi2212 hav- 
ing p ~ 0.2 can be assisted by spin excitations Bear- 
ing this fact in mind, we interpret the peak at 35 mV in 
the difference shown in Fig. 9(a) as the contribution from 
electron tunneling assisted by spin excitations. Indeed, 
the peak position in Fig. 9(a), 35 mV, and the magnitude 
of SG in Fig. 5(b), 2Ase = 34 mcV, have similar values. 

If the interpretation is correct, then the tunneling data 
in Fig. 5(b) present additional evidence that spin fluctua- 
tions mediate the long-range phase coherence in Bi2212. 
First, this shows that charge carriers are strongly coupled 
to spin excitations. Second, by comparing the values of 
the Josephson leRn product of the spectra shown in Figs 
4(b) and 5(b) (7.5 mV and 24.5 mV, respectively), the 
high value of the Josephson product of the spectra in 
Fig. 5(b) can only occur due to the contribution from the 
tunneling assisted by spin excitations because this contri- 
bution is the only difference between the spectra shown 
in Figs 4(b) and 5(b) [see Figs 8(c) and 9(a)]. Conse- 
quently, this signifies that the spin fluctuations mediate 
the phase coherence in Bi2212. 



IX. THE TWO ENERGY SCALES 

The two energy scales shown in Fig. 14, both relate to 
the SC: the linear Ase (or Ap) scale is the in-plane en- 
ergy scale, whereas the scale proportional to Te is the 
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c-axis energy scale, and they have different origins. It is 
important to emphasize that the Age energy scale is ob- 
served by tunneling spectroscopy and ARPES, whereas 
the Ac energy scale is measured in Andreev reflections 
p2[ . However, in some cuprates, the Ac scale is also 
observed in tunneling measurements along the c axis 
p3| , p3|j94| . For example, in near optimally doped YBCO 
single crystals, the maximum magnitudes of the tunnel- 
ing gap into Cu02 planes and along the c axis are 28 meV 
and 19 meV, respectively Q|. The same values measured 
in electron-doped Ndi.85Ceo.i5Cu04_5 (NCCO) under- 
doped in oxygen are 13 meV and 3.5 meV, respectively, 
and the magnitude of the Andreev gap is 3.5 meV [|63[ . 
So, in YBCO and NCCO, tunneling measurements along 
the c axis show the Ac scale. However, it is not the case in 
Bi2212, where tunneling measurements along the c-axis 
show the Asc scale []6|,p5|. 

The symmetries of the two energy scales are discussed 
elsewhere Q: the pairing energy scale, A^c has most 
likely an anisotropic s-wave symmetry like the PG (CDW 
gap) which defines the magnitude of A^c, whereas the 
coherence energy scale, Ac, has the d-wave symmetry. 
Since spin fluctuations mediate the phase coherence, they 
are responsible for the d-wave symmetry of Ac All 
phase-sensitive techniques, obviously, probe the symme- 
try of Ac, and not of Asc, thus they detect the d-wave 
symmetry. At the same time, tunneling measurements 
show a s-wave symmetry of the condensate | |97| , i.e. the 
symmetry of A^c. 

X. THE ANISOTROPY OF THE TUNNELING PG 

Let us concentrate for a while on the tunneling PG. 
The analysis of the data (see above) suggests that the 
tunneling PG is predominantly a charge gap on charge 
stripes in Cu02 planes. This charge gap is most likely 
a CDW gap. It is not difficult to show the shape of 
the CDW gap. As mentioned above, the shape of the 
PG shown in the inset of Fig. 4(b) indicates that the PG 
has most likely an anisotropic s-wave symmetry. Second, 
from the data presented here and in Ref. |32| , the maxi- 
mum magnitudes of in-plane PG and SG in Bi2212 cor- 
relate with each other as Ap^^m ~ 2Asc,m- The intrinsic 
c-axis tunneling data in Bi2212 mesas show that, at low 
temperature, Ap^.c — Asc,c (see Fig.l in Ref. [^). Since 
Asc,c — Asc, in (compare our data with the data in the 
overdoped sample in Ref. Q ) , it is reasonable to assume 
that the value of Ap^^c corresponds to the minimum of 
the in-plane CDW gap because the tunneling along the 
c-axis is preferable from the minimum of in-plane CDW 
gap. Then, one can conclude that the in-plane CDW 
gap in Bi2212 is anisotropic with the anisotropy ratio of 

Aps^raax / Aps jjiiji — 2. 

From quasi- ID structure of the stripes, the CDW gap 
has either two- or four-fold symmetry. The maximum 



magnitude of the SG is most likely defined by the mini- 
mum of the charge gap on the stripes. 

XI. CHANGES IN BI2212 WITH DECREASE OF 
TEMPERATURE 

By synthesizing our findings we briefly describe here 
the changes occurring in Bi2212 with decrease in tem- 
perature. 

Analysis of the data shows that (i) the tunneling char- 
acteristics corresponding to the QP peaks in Bi2212 are 
in excellent agreement with theoretical predictions made 
for a quasi- ID topological-excitation liquid, and (ii) the 
phase coherence along the c-axis in cuprates is estab- 
lished at Tc due to spin fluctuations in local AF domains 
of Cu02 planes. By taking into account (i), we assume 
that the quasi-lD topological excitations reside on quasi- 
ID charge stripes. The assumption is reasonable because 
AF domains which separate the quasi-lD charge stripes 
are truly 2D Consequently, the quasi-lD excitations 
cannot reside into the AF domains. 

Depending on the compound, the charge stripes are 
formed in cuprates at Tcharge 70-300 K (or higher) 
(see references in Ref. |Q ) . The excitations along stripes 
appear with appearance of the stripes, thus, at Tcharge- 
The stripe excitations couple with each other at Tonset < 
Tcharge, which Can be considered as an onset of SC in 
Bi2212. For example, in slightly overdoped Bi2212, two 
independent tunneling studies ||^,^ show that the onset 
of SC occurs at 110-116 K. Finally, the long-range phase 
coherence is established at Tc along the c-axis (i.e. be- 
tween Cu02 planes) due to spin fluctuations. Thus, there 
are 3 characteristic temperatures: Tc < Tonset < Tcharge- 
In addition, between Tonset and Tcharge, there is a char- 
acteristic temperature which was attributed to the 
spin ordering in AF domains . 

The critical temperature, Tc, can be measured directly, 
and Tcharge and Tonset are proportional to Ap^ and Asc 
(see Fig. 7), respectively. It is interesting to note that, 
since Aps ~ 2 Asc, Tonset IS always about twice as smal 
as Tcharge- It is possiblc to estimate the relation between 
Tonset and Asc. By using the data measured in slightly 
overdoped Bi2212 and presented in Refs [0,^, we have 

Tonset ~ 0.4Asc/fc_B. 

Earlier, we proposed a magnetic coupling between 
stripes (MCS) model assuming that the Cooper pairs are 
formed along charge stripes, and the phase coherence is 
established due to magnetic coupling between SC stripes 
pof . From the present work, the general idea of the MCS 
model is correct, however, the more precise name of the 
model of SC in cuprates should be the magnetic coupling 
between striped Cu02 planes. 
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XII. COOPER PAIRS 

As discussed in Section XI, the topological solitons 
condense at Tonset by lowering their energy. Formally, 
the topological excitations shown in Fig. 11 are "pure" 
charge excitations. The question is at what temperature 
is the spin sector condensed, at Tonset or Tc? If Davydov's 
hisoliton scenario |^,^ is realized than the spins con- 
dense together with the charges at Tonset by creating sin- 
glet states. If the condensate consists of kink-kink bound 
states, then spins may condense at Tc- Since "pure" 
charge excitations have boson-like properties, generally 
speaking they do not need to create a bound state. So, 
it is most likely that the Cooper pairs in Bi2212 are 
Davydov's bisolitons, and the phase coherence among the 
bisolitons is established at Tc due to spin fluctuations (ex- 
citations) . 

It is important to emphasize that, in such a scenario, 
formally, the Cooper pairs exist above Tc, between T^ 
and Tonset- However, de fait, they are the local Cooper 
pairs: In order to jump from one stripe to another (above 
Tc), they have similar difficulties as a single electron. 

One would be interested to know what is the glue be- 
tween two electrons in the soliton scenario for the SC in 
Bi2212. It is phonons, and the electron-phonon interac- 
tions are strong and nonlinear |23J23]. 



XIII. DATA OBTAINED IN CUPRATES 

The quasi- ID topological-excitation-liquid scenario is 
attractive, because by using it, it is possible to explain 
experimental data obtained in cuprates. Let's briefly dis- 
cuss some data. 

The discrepancy for the presence of the PC in the over- 
doped region in transport and tunneling measurements 
can be easily understood: they measure two different 
PCs. In tunneling measurements, it is predominantly 
the charge gap, and, in transport measurements, it is the 
resistance to the propagation of stripe excitations along 
stripes which are located in AF environment. For exam- 
ple, if a molecular chain is embedded into a medium, a 
frictional force acts on the soliton . 

In frameworks of the quasi-lD topological-excitation- 
liquid scenario in Bi2212, magnetic and non-magnetic im- 
purities will affect the SC similarly, which is the case in 
cuprates | 100 | 



The oxygen isotope effect [101] can be 
naturally understood since the stripe dynamics above Tc 
depends on the underlying lattice [102|. At the same 
time, a very weak "ECS isotope effect" in cuprates is 
well understood in terms of the soliton SC ||2^,^ . 
Mysterious vortex-like excitations found in the Nernst 



solitons Ig^l ■ Such scenario in Bi2212 can explain why the 
phenomenon of SC in hole- and electron-doped cuprates 
can be understood within a common scheme [^. The 
occurrence of sudden CDW ordering at low temperature 
in undoped cuprates pa^ and the insulating behavior in a 



strong magnetic field |104| correspond to the disappear- 
ance of the stripe excitations. 

By using the quasi-lD topological-excitation-liquid 
scenario in Bi2212 and the fact that the phase coherence 
along the c-axis is established due to spin fluctuations, 
it is possible to explain why Tc and Age do not correlate 
with each other in cuprates: the magnitude of Age is de- 
fined by the magnitude of the charge gap on stripes, at 
the same time, the value of the bulk Tc value depends on 
spin fluctuations into local AF domains. 

It is possible that zero-bias conductance peak observed 
in tunneling measurements l6^ , |9^ co rresponds not only 
to Andreev surface bound states [ 105 1G6| but also to the 
topological solitons [see Fig. 2(a)]. 

The SC condensation energy has the maximum in the 
overdoped region near p = 0.19 |,|ll|. Why? The SC 
in Bi2212 is literally spread into the two channels which 
have different energy scales: in-plane and along the c- 
axis. From Fig. 14, the two energy scales have the same 
values approximately at p = 0.2. Thus, the two energy 
scales emerge into one at p = 0.19-0.2. 

Recently, a new energy scale is found in ARPES spec- 
tra of hole-doped cuprates: an abrupt change (a kink) 
has been observed in the electronic QP dispersion ]107|. 
The kink is associated either with phonons or with the 
magnetic resonance peak. In fact, the clue to the ori- 
gin of the kink in ARPES spectra can be found in Ref. 
]108[. Measurements in 2D layered 2H-TaSe2 CDW com- 



pound foimd a similar energ y sca le in ARPES spectra 
]108]. The authors conclude J 108]: "A reduction in the 
scattering rates below this energy indicates the collapse 
of a major scattering channel with the formation of the 
CDW state accompanying the appearance of a bosonic 
mode in the excitation spectrum of the system." In the 
CDW compound, it is clear that this bosonic mode cor- 
responds to the appearance of solitons in the CDW state. 
Consequently, the origin of the kink in ARPES spectra 
of hole-doped cuprates is the same, i.e. due to solitons. 
In fact, it is also possible that this kink corresponds to 
the real Fermi surface in cuprates and in 2H-TaSe2 (see 
Fig. 2 in Ref. 

Lastly, it is important to note that, in the framework 
of the quasi-lD topological-excitation-liquid scenario in 
Bi2212, one can explain the violation of the sum rule, 
observed in infrared measurements ]109]: the stripe ex- 
citations do not obey the sum rule. 



effect ] |103| correspond to the stripe excitations. The au- 
thors are correct by pointing out that the observed ex- 
citations are vortex-like because vortices themselves are 
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XIV. CONCLUSIONS 
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We would like to emphasize that we do not present 
here an evidence that the tunneling pseudogap in Bi2212 
is a CDW gap. This conclusion is based on the analysis of 
the data. Here we found that the tunneling characteris- 
tics corresponding to quasiparticle peaks are in excellent 
agreement with theoretical predictions made for a quasi- 
ID topological-excitation liquid. In fact, this information 
is sufficient in order to conclude that the tunneling pseu- 
dogap in Bi2212 is predominantly a charge gap. Since 
solitonlike excitations, and not electrons (holes), play 
the role of quasiparticles in Bi2212, this implies that the 
stripes are insulating. In other words, there is a charge 
gap along stripes which is most likely a CDW gap. Then, 
in order to pull out an electron (a hole) from a stripe, it 
is necessary, first of all, to overcome the charge gap on 
the stripe. Consequently, the tunneling pseudogap is pre- 
dominantly the charge gap. However, it is possible that 
the pseudogap shown in Fig. 3(b) consists of two (or a 
few) contributions: the predominant contribution from 
the charge gap and another contribution, for example, 
from AF correlations. 

In summary, tunneling measurements have been car- 
ried out on underdoped, overdoped and Ni-substituted 
Bi2212 single crystals. Tunneling spectra below Tc are 
the combination of incoherent part from the pseudogap 
and coherent quasiparticle peaks. There is a clear cor- 
relation between the maximum magnitude of the pseu- 
dogap and the distance between the quasiparticle peaks. 
The tunneling characteristics corresponding to the quasi- 
particle peaks are in excellent agreement with theoretical 
predictions made for a quasi-lD topological-excitation 
liquid. So, the quasiparticle peaks in Bi2212 appear 
from the condensation of quasi-lD topological excitations 
which reside most likely on charge stripes. Analysis of the 
data suggests that (i) the tunneling pseudogap in Bi2212 
is predominantly a charge gap (CDW gap) on dynami- 
cal charge stripes, which has most likely an anisotropic 
s-wave symmetry, and (ii) the Cooper pairs in Bi2212 
are Davydov's bisolitons. In addition, analysis of data 
measured by different techniques shows that the phase 
coherence is established at Tc due to spin fluctuations 
in local antiferromagnetic domains of Cu02 planes. It 
is possible that the magnon-like excitations which cause 
the appearance of the magnetic resonance peak in inelas- 
tic neutron scattering spectra are in resonance with the 
bound states of solitons |^ . 

Many other data obtained in cuprates can be nat- 
urally understood in the framework of the quasi-lD 
topological-excitation-liquid scenario. In the phase di- 
agram of cuprates, the topological-excitation-liquid con- 
ductor known as a " strange metal" is located between an 
antiferromagnetic Mott insulator at low hole concentra- 
tion and a Fermi-liquid metal at high hole concentration. 
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